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The influence of the -decay radionuclide layer (the energy of x-particles are 5.45 MeV) on the radar cross 
section (RCS) of sphere objects was calculated under different radioactivities, frequencies, and sphere radii. 
When the sphere radius is smaller than 50cm, the tendency of the electron densities of the plasma slab is to 
ascend first and then descend, and the typical maximum electron densities with a radioactivity of 10 Ci/em? 
raises from 7.02 x 10'° to 1.76 x 101! when the sphere radii increases from 10 to 300cm. The average data 
of a normalized RCS of a sphere with radius of 12.5cm, which is coated with a radionuclide layer with dif- 


ferent radioactivities are —0.35, —0.50, 


respectively. 


0.79 and —1.13 dB when the radioactivity is 1, 2, 5 and 10 Ci/em’, 
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I. INTRODUCTION 


In the recent years, the plasma stealth technique has gotten 
more and more attentions because of its smaller impact on 
the shape of objects and better suitability in a wide frequency 
band [1-3]. Many techniques can produce plasma: arc dis- 
charge, laser excitation, electron guns, or radioactive iso- 
topes [4-7]. When the frequencies of plasma generated from 
the radioactive isotopes commensurate with the frequencies 
of electromagnetic waves (EMWs), the damping loss result- 
ing from the excitation in the plasma slab by EMW reaches 
the maximum. Therefore, the plasma stealth technique by 
radioactive isotopes is superior, especially in the ultra high 
frequency (UHF) and super high frequency (SHF) bands. 


Sphere metal objects are the most useful in the calibration 
body of the radar cross section (RCS) field [8, 9]. In past 
years, many teams had worked out the influence of RCS on 
different plasma densities and shapes, but most of their re- 
search was under their own definition of density or shape [10- 
17]. These works could be considered a verification of plasma 
theory, but for actual application, there is still a pronounced 
gap. This paper aims to calculate the influence of plasma in- 
duced by a radionuclide layer (RNL), which is coated on the 
surface of a sphere metal object, on the RCS in different ra- 
dioactivities and ratios of the radii of the sphere to the elec- 
tromagnetic wavelength (RRW). 
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Il. THEORETICAL PRINCIPLE AND RESEARCH 
METHODS 


A. Model 


Figure 1 is a schematic drawing of the simulation. A 
thin RNL is coated on the surface of the spherical metal ob- 
ject. The common radioactive isotopes of «-decay are Po-210 
(5.304 MeV, 138.4d), Pu-238 (5.456 and 5.499 MeV, 87.7a) 
and Am-241 (5.443 and 5.485 MeV, 432.6a). The energy of 
the -particles ejected from all the three isotopes are similar, 
meaning the energy loss of a distance unit of these «-particles 
is also similar. In our paper, all the calculations are under the 
x-particle energy of 5.45 MeV as the representative, which 
possesses no obvious differences of energy loss in air with 
the three aforementioned isotopes. We assume the material 
of the spherical metal object is the prefect electric conductor 
(PEC) because the surface effect depths of most metals are 
much less than the regular radius of the sphere object. When 
the EMW propagates into the plasma slab referable to the 
ionization effect of «-particles, the damping loss of charged 
particles would decrease the energy of the EMW, and it also 
would refract the spatial distribution of EMW. All these will 
influence the RCS of the spherical metal objects. 
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Fig. 1. (Color online) The schematic drawing of the simulation. 
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B. Electron density distribution in plasma slab induced by 
x-particles 


As the velocity of x-particles is much less than light, we 
can get the average number of new-born ions and ioniza- 
tion electrons by a single «-particle using the nonrelativis- 
tic Bethe-Bloch formula and an average ionization energy of 
36.08 eV in air [18] 


— 204 2 
( =) X 4nz*e NE m (4 | (1) 


dx (mov?) I 


The frequency of plasma is mainly determined by the den- 
sity of electrons, because the mass of the electron is much 
less than the ion. In the simulation, we could only trace these 
ionization electrons, as a result. In the air, the new-born ion- 
ization electrons mainly have three physical processes: dif- 
fusion, recombination with positive ions, and adsorption by 
neutral molecules. During these processes, the negative oxy- 
gen ions formed by adsorption dissociate rapidly under ul- 
traviolet rays when the ionized electrons reach a specific en- 
ergy [19]. In our calculation, we neglected the adsorption 
process and the elimination would not bring any substantive 
effects to the results. 

The diffusion velocity of electrons is related to the gradi- 
ent of the electron density (ne), and considering the effect of 
ambipolar diffusion in the plasma slab, the coefficient of dif- 
fusion is Dait = 3.7 x 10? cm?/s. 


dN 
—) = Dä . 2 
( dt i dit V Ne ( ) 


The recombination velocity is proportional to the concen- 
trations of both positive and negative ions. The whole plasma 
slab system is close to neutral, so we could consider that the 
numbers of positive ions and electrons are nearly the same, 
and the coefficient of recombination is a = 1.5 x 107° cm/s 
(both the coefficients are validated by experiments). 


(=) : =ant ne & anz. (3) 


The final density of electrons is determined by diffusion, 
recombination, and ionization by «-particles. Figure 2 is the 
flow chart of the simulation. The time gap of each cycle is the 
average time gap between the adjacent two ejected «-particles 
in one square centimeter, which is determined by the radioac- 
tivity of the RNL. 


C. The interaction between EMW and plasma 


In the plasma slab, there exist four kinds of collisions be- 
tween electrons and molecules, electrons and ions, ions and 
molecules, and electrons and electrons. The collisions mostly 
happen between electrons and molecules in the troposphere. 
Ginzburg defined a coefficient named effective collision fre- 
quency [20], when the frequency of EMW is much smaller 
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Fig. 2. (Color online) The flow chart of the simulation of electron 
density. 


than the effective collision frequency, the latter could be for- 


mulated as 
Nm IT 
— 4| —_ 4 
2.7 x 1019 V 300° @) 


where Nm is the number of molecules per unit of volume. 

But the effective collision frequency exhibits an impact on 
the EMW frequency. The relationship with the relative di- 
electric constant and effective collision frequency is 
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where Kı and Ko are two coefficients related to ratio of 
EMW frequency to the effective collision frequency [20], and 
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The plasma slab induced by the RNL, which is coated on 
the surface of the spherical object, is isotropic, that means the 
density distribution is only related to the radial distance, r. So 
the plasma slab could be divided into m layers and the radius 
of the boundary between the i — 1" and i layer is set as aj: 
when 0 < r < aj, the material is a prefect electric conductor 
(PEC); the relative dielectric constant of the i layer is Eri 
and the relative permeability is 1, approximately. The rela- 
tionship between the RCS and scattering coefficients b» and 
Cn 1S 
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The scattering coefficients bn and cn are calculated by 
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And the recurrent relationship of the coefficient Ry, and 
Re m 
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And the monostatic RCS is one of the particular cases (0 = 
7), so the RCS could be simplified as 


2 [0.6] 
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n=1 


All the RCS values are represented in the normalized re- 
sults, 101g as in which r is the radius of the spherical 
object. 
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Il. RESULTS AND DISCUSSION 
A. Electron density distribution in plasma slab 


The average range of 5.45MeV «-particles is 3.84cm in 
the air at normal pressure and temperature; the track of a sin- 
gle «x-particle is almost straight along the ejection direction. 
By the simulated average energy loss along the radial direc- 
tion in consideration of the surface curvature of the spherical 
metal objects and the randomness of ejection direction, we 
can get the average electron density per unit volume ionized 
by a single «-particle, shown in Fig. 3. The solid angle per 
square centimeter on the surface increases with the decrease 
of sphere radius, so the average electron densities ionized by a 
single x-particle decrease. When the sphere radius is smaller 
than 50 cm, the results show that the tendency of electron den- 
sities is to ascend first and then descend. It could be attributed 
to the effect of the sphere curvature for a relatively small ra- 
dius, which would cause the decrease of the effective RNL 
surface with distance, in the proximal region to the surface. 
Under the combined interaction of the effective RNL surface 
and the volume corresponding to a unit of solid angle, the ra- 
dial electron densities possess tendency phenomena with the 
distance to the sphere. 
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Fig. 3. (Color online) The average radial electron density by a single 
x-particle per square centimeter. 


Based on the average radial electron density and flow chart 
shown in Fig. 2, we have calculated the actual electron den- 
sity at different radii and radioactivities. In the calculation, 
we set the thickness of the plasma slab at 50cm, which is 
much larger than the range of «-particle. The consideration 
of the diffusion effects and the results are shown in Fig. 4. 
The electron density increases with the radius when the ra- 
dioactivity of RNL is 1 uCi/cm?; the maximum density data 
are 2.42 x 107, 3.52 x 107, 4.48 x 107 and 5.17 x 107 cor- 
responding to the sphere radii of 10, 25, 50 and 100 cm, re- 
spectively, and the results are close to a value 5.76 x 107 for 
an infinite large plate (r — 00). 

When the radioactivity increases 10 times, the ratios of the 
electron density data with radioactivity at a certain value to 
the data with ten times the amount of radioactivity are ba- 
sically consistent for the same sphere radii. The ratio, ex- 
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Fig. 4. (Color online) (a) The electron density of 1 uCi/cm? with distance to sphere surface at different radii, (b) the activity coefficients at 


different radii. 


pressed as a(x) and denominated by the activity coefficient, 
is shown in Fig. 4(b) at different sphere radii. Then we sum- 
marized an empirical formula to describe the electron density 
in the plasma slab induced by 5.45 MeV «-particles with the 
distance to the sphere surface. When the radioactivity is in 
the range of 1 uCi/cm? to 10 Ci/cm?, the differences between 
the empirical equation and the simulation results are less than 
3%. 


S (15) 


n(x) = no(x)a(a) 
where no(x) represents the electron density data at 1 wCi/cm? 
in Fig. 4(a), x is the distance to sphere surface, A is the sur- 
face radioactivity in unit of uCi/cm?. 
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Fig. 5. (Color online) The electron density with x at 10 uCi/cm? of 
radioactivity and different radii by the empirical formula. 


In accordance with the empirical formula, we can get the 
distribution of electron density in conditions of different radii 
and radioactivities conveniently. Figure 5 is the electron den- 
sities with different sphere radii when the radioactivity is 
10 Ci/cm? by the empirical formula. When the sphere radii 
are small, the regions of maximum electron density are not 
close to the sphere, the tendency of electron density is ascend 
to first and then descend. For example, when the sphere ra- 


dius is 10cm and the radioactivity of RNL is 10 Ci/cm?, the 
distance between the region of the maximum electron density 
and the sphere surface is about 1 cm. 
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Fig. 6. (Color online) The monostatic RCS at different radioactivi- 
ties. 


B. The influence on monostatic RCS of «-decay RNL 


One of the most important parameter is the RRW in the 
monostatic RCS research on spherical objects. Normally, the 
normalized RCS data is completely consistent for the same 
RRWs, due to the relativity of the spatial scale. But the dis- 
tributions of electron density in the plasma slab outside the 
spherical object will change with the radii and radioactivity. 
So in our model, even if the RRWs are the same, there ex- 
ist differences in the normalized RCS data. Figure 6 is the 
monostatic RCS data in the radioactivity range from 1 Ci/cm? 
to 10Ci/cm? and at a radius of 12.5cm. The amplitudes 
of oscillation decrease with the increase of RRW and when 
r/Xo > 1, the average monostatic RCS data of spherical ob- 
jects without RNL is 0. The average values decrease with 
the increase of the radioactivity of RNL: the average data is 
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—0.35, —0.50, —0.79 and —1.13dB, respectively, when the 
radioactivity is 1, 2, 5 and 10 Ci/em?. Also, the amplitudes 
of oscillation would decrease with the increase of radioactiv- 
ity, and there exists deviations in the positions of peaks and 
valleys. The reason is the refraction of incident and scatting 
waves caused by the inhomogeneous distribution of plasma. 
The higher electron density indicated by higher radioactivity 
would result in more obvious refraction. 
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Fig. 7. (Color online) The monostatic RCS at different sphere radii. 


Figure 7 shows the monostatic RCS data at different sphere 
radii and a radioactivity of 10 Ci/cm?, and we can also find 
the differences in these curves with different sphere radii for 
the same radioactivity. Similarly, the amplitudes of oscilla- 
tion decrease with the increase of RRW or the decrease of the 
sphere radius. Certain deviations also exist in the positions of 
peaks and valleys. The average values of normalized mono- 
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